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ABSTRACT 

O , The global temperature distribution of the cosmic gas-phase oxygen at z—3 is de- 

termined by combining high resolution cosmological simulations of individual proto- 
galactic, as well as larger, regions with the observed, extinction-corrected, rest-frame 
V-band galaxy luminosity function. The simulations have been performed with three 
different stellar initial mass functions (IMFs), a Kroupa (K98), a Salpeter (S) and an 
\ Arimoto- Yoshii (AY), spanning a range of a factor of five in chemical yield and spe- 

J> ■ cific SNII energy feedback. Gas-phase oxygen is binned according to T as log(T) <4.0 

OO i ("cold"), log(T) - 4.5 ("warm"), and log(T) -5.0, 5.5, 6.0, 6.5, 7.0 ("hot" phases). 

G\ ■ Oxy gen is found to be distributed over all T phases, in particular for the top-heavy AY 

IMF. But, at variance with previous works, it is found that for the K98 and S IMFs the 
cold phase is the most important. For these IMFs it contains 47 and 37%, respectively, 
of all gas-phase oxygen, mainly at fairly high density, nn> 0.1 cm -3 . The implications 
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, of this in relation to observational damped Lyman-a absorber (DLA) studies are dis- 

cussed. In relation to "missing metals" it is found that a significant fraction of the 
oxygen is located in a warm/hot phase that may be very difficult to detect. Moreover, 
it is found that less than about 20-25% of the cosmic oxygen is associated with galax- 
Q ■ ies brighter than My — —22, i.e., the faintest galaxy luminosities probed by current 

\ metallicity determinations for Lyman Break Galaxies (LBGs). Hence, 75-80% of the 

oxygen is also in this sense "missing". From the LBG based, A— 1500 A UV luminosity 
density history at z>3, we obtain an essentially IMF independent constraint on the 
mean oxygen density at z=3. We compare this to what is obtained from our models, 
for the three different IMFs. We find that the K98 IMF is strongly excluded, as the 
chemical yield is simply too small, the Salpeter is marginally excluded, and the AY 
matches the constraint well. The K98 IMF can only match the data if the A— 1500 A 
extinction corrections have been overestimated by factor of —4, which seems highly 
unlikely. The yields for K98 are also far too small to match the observational data for 
C iv. The optimal IMF should have a yield intermediate between the S and AY. 
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1 INTRODUCTION 

A number of "problems" have been discussed in relation 
to the cosmic mean metallicity and/or mean metal density. 
Pagel (1999) discussed the low redshift (z~0) "excess metals 
problem" showing that a comparison of the average metal to 
stellar densities in the low- z Universe indicates that the aver- 
age cosmic stellar initial mass function (IMF) has a higher 
yield, than what is obtained for the "standard" Salpeter 
(1955) IMF. This is a well known result for ellipticals in 
clusters (e.g., Romeo et al. 2006 and references therein), but 
Pagel's analysis indicates that it could be more universal. 

Conversely, Pettini (1999) formulated the high-z "miss- 



ing metals problem" as follows: Studies of the co-moving 
rest-frame UV luminosity density of high-z Lyman break 
galaxies (LBGs) allow us to trace the cosmic star formation 
density (or history, SFH), p\(z), up to redshifts z max w 7—8. 
Assuming an IMF of such stars, one can compute the specific 
fraction of heavy elements ("metals") they produce, y, and 
derive the metal production rate pz(z) = yp+(z). The inte- 
gral from z m ax gives, at any given z, the density of cosmic 
metals p% F (z), or, expressed in units of the critical density, 
Q Z FH . Moreover, if one restricts the analysis to elements, 
such as the a-elements, produced almost exclusively in mas- 
sive stars undergoing core collapse, the apparent dependence 
of Pz FH ( z ) °n the IMF is essentially removed, since both the 
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UV light and the oxygen production originate from mas- 
sive stars. This will be discussed further in Section 2, and 
was already noticed by, e.g., Songaila et al. (1990), Madau 
et al. (1996), Pagel (1999) and Pettini (1999). 

Early searches in cosmic structures for which the 
metal/baryon mass ratio (metallicity, Z = Q,z/Qb) can be 
derived either via intergalactic gas quasar absorption line 
experiments (DLAs or the Lya "forest") or through direct 
spectroscopic studies of Lyman break galaxies have found 
that only Sl° BS < 0.20fi§ FH is stored in these components, 
i.e., the large majority of the metals are "missing". Simi- 
lar missing metal problems have been formulated by Wolfe 
et al. (2003) and Prochaska et al. (2003, 2006) on the basis 
of star formation rates and metallicities of DLAs alone. 

Ferrara et al. (2005, FSB05) attempted to quantify 
more precisely the extent of the missing metals deficit, 
and also to suggest where the "missing metals" might be 
found. They found from considering typical stellar masses 
and metallicities of Lyman break galaxies, and an esti- 
mate of the co-moving density of LBGs a contribution of 
qLbg _ g ^ x f rom L yman break galaxies. They also 

estimated the contribution from metals in DLAs to about 
Q.z LA = 3.8 x 10~ 7 . Using the SFH estimate of Bouwens 
et al. (2004b), assuming a Salpeter IMF, and applying a 
dust correction factor of 4.5 (Reddy & Steidel 2004) they 
estimated that the Universe should be characterised by a 
total metallicity density of Q.% FH = 1.84 ± 0.34 x 10~ 5 at 
2=2.3. From the above numbers they concluded that about 
80% of the metals in the Universe are missing, in the sense 
that they are not directly associated with the above two 
galactic components. 

FSB05 suggested that the "missing metals" are associ- 
ated with galaxies, mainly residing in their "hot" gas ha- 
los, having been deposited there by star-burst driven super- 
winds (see also Pettini 2004). This hypothesis seems rea- 
sonable, given that outflow velocities of 300-400 km/s are 
routinely inferred from the spectra of LBGs (e.g., Pettini 
et al. 2001, Shapley et al. 2003). Moreover, semi-analytical 
as well as fully hydrodynamical models of galaxy formation, 
based on Cold Dark Matter (CDM), require such "feed- 
back" in order to enable the formation of realistic galax- 
ies, solving the "over-cooling" , "angular momentum" , "miss- 
ing satellites" and other problems (e.g., Sommer-Larsen 
et al. 1999, Cole et al. 2001, Thacker & Couchman 2001, 
Sommer-Larsen et al. 2003, SGP03). 

Considering the widely used probes of Lya forest metal 
absorption in QSO spectra, Civ and Ovi (e.g., Songaila 
2001, Bergeron et al. 2002, Carswell et al. 2002, Boksen- 
berg, Sargent & Rauch 2003, Schaye et al. 2003, Simcoe 
et al. 2004, Aracil et al. 2004, Scannapieco et al. 2005; and 
Bergeron & Herbert- Fort 2005, Tripp et al. 2006), FSB05 
inferred observed integrated C IV and O VI column densities 
over the redshift range 1.7<2<3.8, and the corresponding 
cosmic average densities of these. Assuming a two-phase 
intergalactic medium (IGM), consisting of a "cold" part 
(T~10 4 K), responsible for the hydrogen Lya forest absorp- 
tion, and a "hot" component (T~10 6 K), they showed, using 
a number of simplifying assumptions, that more than 90% 
of the metals can reside in the hot component, without vi- 
olating the C IV and O VI Lya forest metallicity constraints 
(and with the cold phase as the major contributor to the 
C IV and O VI column densities). This result follows provided 



that That ~10 5,5 K, and the density of the hot metal-enriched 
gas is only factors of a few times the mean cosmic bary- 
onic density, which, e.g., at 2=3, corresponds to a hydrogen 
number density of 71h~1.2x 10~ 5 cm~ 3 . FSB05 suggest this 
metal containing gas to be identified as wind-blown galaxy 
halo gas. 

The low luminosity galaxies at high redshift are traced 
by the DLAs (e.g. Fynbo et al. 1999; Haehnelt et al. 2000, 
Schaye 2001a; M0ller et al. 2002). From very early on it was 
found that DLAs contain different regions with very dif- 
ferent temperatures and ionisation states (Turnshek et al. 
1989). The DLA metal content discussed by FSB05 only 
accounts for the cold, mainly neutral component. The Civ 
and Si IV lines in DLAs corresponds to a warmer phase with 
more turbulent kinematics. The cross-section for strong C iv 
and Si iv absorption is much larger than for DLAs consis- 
tent with the picture that this gas is located in an extended 
wind-blown halo around the stellar components (Petitjean 
& Bergeron 1994, see also Adelberger et al. 2003, 2005, Por- 
ciani & Madau 2005, Scannapieco 2005). An even hotter 
component, traced by Ovi and Nv (temperatures of or- 
der 3 x 10 5 K), has recently been identified in DLAs (Fox 
et al. 2007). Fox et al. find that if the temperature of the 
O vi bearing gas is 10 6 K or higher, then this hot phase can 
contribute significantly to the metal mass budget — quali- 
tatively consistent with the proposal of FSB05. 

Bouche et al. (2006a,b) discussed the missing metals 
problem at 2^2.2-2.5 and found that the problem is not 
quite as severe as at 2~3. Counting metal contributions 
from stars in BX galaxies (the equivalent of LBGs at 2=2.2; 
e.g., Adelberger et al. 2004) and DRGs ("distant red galax- 
ies"; e.g., Franx et al. 2003, van Dokkum et al. 2003), as 
well as from gas in SMGs (sub-millimeter galaxies; e.g., 
Blain et al. 2004, Greve et al. 2005) and DLAs (e.g., Pet- 
tini et al. 2003) they could account for about 1/3 of the 
metals expected. The BX galaxies contain the largest frac- 
tion of the identified metals, about 55%. Moreover, Bouche 
et al. (2007) found that another about 1/3 of the metals ex- 
pected can be identified in the IGM at such redshifts. 

Recently, Dave & Oppenheimer (2006, DO06) mod- 
eled the enrichment history of the Universe and address 
the missing metals problem using a fully numerical ap- 
proach. They perform a moderate resolution cosmological 
hydro/gravity simulation of a cubic region of the Universe 
of 32 /i _1 Mpc box size. Their simulation invokes metallic- 
ity and UV background dependent radiative cooling, star- 
formation and chemical evolution, in the instantaneous re- 
cycling approximation, and based on the Salpeter IMF. 
The simulation resolves galaxies of stellar masses down to 
~10 9 Mq, corresponding to a V-band absolute magnitude 
My~ —20 at z~3. Due to the limited resolution, the authors 
adopt a parameterized description of star-burst driven out- 
flows in the form of a "momentum-driven" wind, found by 
Oppenheimer & Dave (2006, OD06) to yield the best match 
of the results of their simulations to various observational 
data. 

DO06 quantify their results in terms of cosmic metal 
fractions in five phases: a) stars, b) star-forming gas, c) 
halo gas (gas inside of the virial radius of galaxy halos, 
which is not star- forming) , d) shocked IGM (gas outside of 
galaxy halos of T>3xl0 4 K), and e) diffuse IGM (gas out- 
side of galaxy halos of T<3x 10 4 K). At z~3, they find that 
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the dominant metal phase is the diffuse IGM, which con- 
tains about 40% of the cosmic metals, with the other four 
phases containing approximately equal fractions of about 
15%. Hence, DO06 find that indeed only about 20% of the 
metals in the 2~3 Universe reside in stars, while the "miss- 
ing" about 80% are located in the gas phase. Moreover, they 
find that the hot halo gas is not the dominant, metal con- 
taining phase. Instead, a large fraction of the missing metals 
are "hidden" in the diffuse IGM. DO06 suggest that this is 
possible, in relation to observations of, e.g., low density IGM 
C IV abundances, because the diffuse gas is somewhat hot- 
ter than assumed in previous estimates, which at gas over- 
densities of 8 = pgas/pgas~ 1 — 100 implies larger Civ to C 
ionisation corrections, and hence that larger amounts of C 
can be "hidden" in this phase. 

In this paper we take another approach from that of 
DO06 (see also Calura & Matteucci 2004, 2006). As we will 
show in the paper, due to resolution limitations DO06 likely 
fail to account for about half of the metals in the Universe 
produced by z~3, on top of which adds cosmic variance ef- 
fects, given the relatively small computational box of DO06. 
To undertake simulations of larger cosmological volumes, at 
the same time probing 8-10 magnitudes deeper (Sec. 5) is at 
present, as well as in any foreseeable future, computationally 
prohibitive. Using K-band observations of a sample of LBGs, 
Shapley et al. (2001) determined the rest-frame V-band 
galaxy luminosity function. In Sommer-Larsen & Fynbo 
(2007, paper I) we correct the Shapley et al. luminosity func- 
tion for extinction and other effects, to obtain a "true" z~3 
Lyman break galaxy LF. Subsequently, cosmological high- 
resolution hydro/gravity simulations of the formation and 
evolution of individual galaxies are combined with the cor- 
rected galaxy luminosity function to obtain estimates of the 
average cosmic density of metals (in particular oxygen) re- 
siding in stars at z=3. We consider models based on three 
different stellar IMFs, spanning almost a factor of five in 
chemical yield as well as thermal/kinetic energy feedback 
from supernova type II (SNII) explosions. These are the 
Kroupa (1998, K98) IMF, a typical IMF suited for chem- 
ical evolution models of the Solar Neighbourhoohd, e.g., 
Boissier & Prantzos (1999), the "standard" Salpeter (1955, 
S) IMF and the Arimoto-Yoshii (1987, AY) IMF, which is 
well suited for describing the chemical evolution of elliptical 
galaxies. The simulations have sufficiently high resolution to 
allow a two-phase modeling of the star-forming interstellar 
medium (ISM), consisting of a "cold" T~10 4 K star- forming 
phase and a "hot" T~10 5 - 10 6 K phase, intermixed with the 
cold gas. Star-bursts drive galactic winds by depositing ther- 
mal energy from multiple SNII explosions in the gas, part 
of which is subsequently converted into kinetic energy self- 
consistently by the hydro-code. 

In paper I we conclude that for none of the IMFs is it 
possible to reconcile the amount of oxygen locked in stars 
with the amount predicted from the (observed) cosmic UV 
luminosity density history, hence confirming the "missing 
metals problem" as stated by Pettini (1999) and FSB05. 

In this paper we combine the high-resolution galaxy for- 
mation models with the V-band luminosity function to de- 
termine the cosmologically averaged amount and properties 
of z~3 gas-phase metals, mostly focusing on oxygen, but also 
on carbon, in particular in relation to QSO absorption line 
determinations of the cosmic C IV density, as well as iron, 



in relation to DLA abundances. Combining the results ob- 
tained with those of paper I the total cosmic metal distribu- 
tion is obtained for the three IMFs considered. Comparing 
in turn these results to inferences from the observed cos- 
mic UV luminosity density history allows us to significantly 
constrain the "true" cosmic stellar IMF. Finally, relations 
to the "missing metals problem" , as well as results obtained 
by other authors, are discussed. 

The paper is organised as follows: In section 2 we derive 
constraints on cosmic metal production from the cosmic UV 
luminosity density history, in section 3 we briefly describe 
the hydro/gravity galaxy formation simulations, and in sec- 
tions 4 and 5 we present the approach used in this paper to 
determine the temperature distribution of the cosmic gas- 
phase metals. Section 6 presents results on gas-phase C IV 
abundances, and, relating these to observations of C IV ab- 
sorption lines in QSO spectra, we derive a constraint on the 
cosmic IMF. In section 7 we combine the results of paper 
I and those obtained here to derive the total cosmic metal 
distribution. Comparing this to what is obtained from the 
cosmic UV luminosity density history we obtain an addi- 
tional constraint on the cosmic stellar IMF. In section 8 wc 
relate our results to the "missing metals problem" and com- 
pare them to those of other workers in the field, in section 9 
we demonstrate that our results are robust to changes of the 
numerical resolution and, finally, section 10 summarises our 
conclusions. 

In the paper we assume the flat A cosmology, with Qm 
= 0.3 and S1 A =0.7, and H = WOh km/s/Mpc, with ft=0.7, 
unless it is explicitly stated otherwise. 



2 CONSTRAINTS ON METAL PRODUCTION 
FROM THE COSMIC UV LUMINOSITY 
DENSITY HISTORY 

Constraints on metal production from the cosmic UV lu- 
minosity density history are discussed in paper I, but for 
coherence we present the main discussion in the following as 
well. 

The luminosity of young galaxies at rest-frame A ~1500 
A is a measure of the rate of formation of massive 
stars (mainly O and B type) as shown by, e.g., Madau 
et al. (1998). Moreover, a-elements like oxygen are almost 
exclusively produced in such massive stars, and hence there 
is a direct link between the cosmic average oxygen pro- 
duction rate density and the average UV luminosity den- 
sity, which is essentially independent of the stellar ini- 
tial mass function — see, e.g., Pettini (1999). To estimate 
the cosmic oxygen production rate density as a function 
of redshift we use various recent estimates of the average 
cosmic star formation rate. The estimates have been ob- 
tained from Steidel et al. (1999), Bouwens et al. (2004a,b), 
Bunker et al. (2004) , Giavalisco et al. (2004) , Schimonovich 
et al. (2005), Bouwens et al. (2006) and Sawicki & Thomp- 
son (2006). The estimates used are all based on assuming a 
"standard" Salpeter IMF in converting from UV luminosity 
to star formation rate. The oxygen (mass) yield of a Salpeter 
IMF is 0.01 (e.g., Lia et al. 2002a,b) — hence in Fig. 1 we 
have multiplied the published SFRs by 0.01 to obtain the 
(essentially) IMF independent oxygen production rate den- 
sities. The estimates have moreover been corrected to cor- 
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Figure 1. Cosmic oxygen enrichment history based on ob- 
servations of high-2 Lyman break galaxies. The data points 
shown are based on Steidel et al. (1999; green crosses), Schimi- 
novich et al. (2005; open red dot), Sawicki & Thompson 
(2006; open blue dots), Giavalisco et al. (2004; black squares), 
Bunker et al. (2004; blue square), Bouwens et al. (2004a; ma- 
genta triangle), Bouwens et al. (2006; magenta circle), Bouwens 
et al. (2004b; magenta pentagon — note that the latter only rep- 
resents L >0.3xL^j V z _ 3 ). The "maximum" model is shown by 
the upper short-dashed curve, the "median" model by the long- 
dashed, and the "minimum" model by the lower short-dashed. 

respond to (at any redshift) the UV luminosity density of 
galaxies brighter than 0.1xL^ z=3 following an approach 
similar to that of Bouwens et al. (2006; note though that 
the 2~7.5 data point of Bouwens et al. 2004b only repre- 
sents L >0.3xLy Vz=3 ). Finally, the values shown in the 
plot result from multiplying the observed values by a dust- 
correction factor of 5.5. This value is intermediate between 
the z ~3 values of 4.5 and 6.5 suggested by Reddy & Stei- 
del (2004) and Dahlen et al. (2007), respectively (but see 
below) . 

The cosmic average oxygen density at 2=3 can now be 
obtained by integrating the oxygen production rate density, 
po(z), from 2=3 and back in time. We shall consider three 
models for po{z): the "median" model, the "maximum" 
model, and the "minimum" model. The median model is ob- 
tained by calculating the median values of po in the bins 2=3 
to 4 and 2= 4 to 6 (assigning equal weight to each data point, 
but excluding the 2=5.9 value of Bunker et al. 2004; see be- 
low), and connecting these values with the z ~7.5 value of 
Bouwens et al. (2004b). This model is shown in Fig. 1 by 
the long-dashed curve. The maximum model is obtained by 
connecting the 2=3 to 4 median value with the Giavalisco 
et al. (2004) z ~5.7 value, and the Bouwens et al. (2004b) 
z ~6.0 value, and multiplying the resulting "upper enve- 
lope" by a factor 6.5/5.5 to maximize also the dust cor- 
rection. This model is shown by the upper short-dashed 
curve in Fig. 1. Finally, the minimum model is obtained by 
connecting the 2=3 to 4 median value with the Giavalisco 
et al. (2004) z ~4.8 value, the Bouwens et al. (2004a) value 
and the Bouwens et al. (2006) z ~7.5 value, and multiply- 



Figure 2. Relation between virial radius and absolute V-band 
(rest-frame) magnitude at 2=3 for the K98 (green), Salpeter 
(blue) and AY (red) IMFs. 

ing this "lower envelope" by a factor 4.5/5.5 to minimize the 
dust correction (but see below) . This model is shown by the 
lower short-dashed curve in Fig. 1. 

Assuming a flat space world model, the 2=3 cosmic av- 
erage oxygen density can now be evaluated as 

po(2=3) = - f po(z(t)) dt = 

J t(z— 3) 



H ° J 3 (1 + z)y/(l + 2) 2 (1 + zQ M ) ~ 2(2 + z)Q A 

where values of Ho = lOO/i km/s/Mpc, with h=0.7, Qm = 
0.3 and Qa~0.7 are assumed in this paper. Using eq. 1 we 
obtain po(z=3) = 0.32, 0.42 and 0.56 xlO 7 Af Q /(/i- 1 Mpc) 3 
for the "minimum" , "medium" and "maximum" models, re- 
spectively. These numbers represent an integral constraint 
that any successful model of galaxy formation must meet. 
As will be shown in section 5 this enables us to strongly 
constrain the "cosmic" stellar IMF at z >3. 

In determining the "median" and "minimum" models 
we excluded the 2=5.9 value of Bunker et al. (2004). This is 
simply done on the basis of the large discrepancy between 
this value, and all other 2=4 to 6 values. To indicate the 
effect of this, we determined an alternative median model 
in which the Bunker et al. (2004) value is included in the 
2=4 to 6 bin. This changed the median model estimate of 
po(2=3) from 0.42 to 0.40 x 10 7 Me/^^Mpc) 3 , i.e. a ~2 % 
change, and hence quite small effect, which, given all other 
uncertainties, we shall ignore in the following. 

Bouwens et al. (2006) propose that the A -1350-1500 A 
extinction correction decreases with increasing z from z ~3. 
Converting their proposed, extinction-corrected SFR(2) to 
a corresponding po (2) (including changing the UV luminos- 
ity density from their limit of 0.04xLy Vz=3 to the limit 
of 0.1xL£, Vz=3 used here) results in a value of po{z—3) = 
0.33 x 10 7 M©/(/i _1 Mpc) 3 , hence within the bounds derived 
above. 
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Figure 3. Cumulative mass of gas-phase oxygen with galacto-centric distance for the same three galaxies at z=3, simulated with the 
K98 IMF (left) and the AY IMF (right). Gas temperature colour coding is as follows: log(T)=4: green, 4.5: blue, 5: cyan, 5.5: magenta, 
6: black, 6.5: red, 7 black (thick line). For all galaxies the region from r v i r to 8 r v i T is shown. For the K98 IMF, the galaxies have 
V-band absolute magnitudes of —16.1, —21.4 and —25.3. For the AY IMF, the corresponding numbers arc —13.8, —19.3 and —25.4. 
The corresponding stellar masses are 5.8xl0 7 M , 7.3xl0 9 A/ Q and 3.4X10 11 AIq for the K98, and 6.4xl0 6 M©, 5.2x10 s M and 
2.5X10 11 Mq for the AY IMF. 



The integral constraint obtained in this section can 
also be expressed in units of the critical density: we obtain 
n (z=3) = po(z=3)/p clit = 0.81, 1.05 and 1.41 xlO" 5 for 
the minimum, median and maximum models, respectively. 

FSB05 obtained for the total density of heavy elements 
in units of the critical Q,z = (1.84±0.34) x 10~ 5 , by integrat- 
ing the observed cosmic star formation history (SFH, based 
on a Salpeter IMF), as presented by Bouwens et al. (2004b), 
to .z=2.3 and applying a dust correction factor of 4.5. The 
above value translates into fi (>=2.3) = (0.77±0.14)xl0~ 5 . 
If the oxygen density production rate models shown in Fig. 1 
are continued 2=2.3, we would obtain 1.22, 1.55 and 2.01 
xlO" 5 . The reason for the apparent discrepancy is two- fold: 
a) we use a dust correction of a factor of 5.5, and b) more 
importantly, the Bouwens et al. (2004b) SFH was based on 
a UV luminosity to a limit of 0.3xL^ Vz=3 , rather than the 
limit of O.lxLuy z=3 adopted in this work. Using the z ~ 3 
UV luminosity function of Adelberger & Steidel (2000) we 
estimate the latter correction to be a factor of 1.5. Mul- 
tiplying this by a factor of 5.5/4.5 the findings of FSB05 
would translate into about (1.4±0.3) x 10 -5 , in good agree- 
ment with our values above. 

Extremely dust-obscured galaxies such as high-z 
SCUBA sources (e.g., Smail et al. 1997; Barger et al. 1998; 
Eales et al. 1999; Chapman et al. 2005) are not included in 
the dust-corrected star formation rate densities shown in 
Fig. 1. Most such galaxies have UV luminosities that heav- 
ily under-represent their star formation rates and thus such 
galaxies are not properly included in the estimate of the 
oxygen production rate density. Their numbers may be suf- 
ficiently large that their contribution is significant, however 
we neglect this contribution for two reasons: a) the contri- 
bution from such galaxies would be largest at z~2, where 
the redshift distribution of SCUBA sources appears to peak 
(Chapman et al. 2005, Reddy et al. 2005), whereas we are 



concerned with the range 2 > 3, and b) more importantly, 
we are concerned with the oxygen production rate density 
history of Lyman Break Galaxies only, as we wish to com- 
pare the integral constraint thus obtained to results of com- 
bining the observed LBG 2~3 optical rest-frame luminosity 
function with detailed, high-resolution models of galaxy for- 
mation. 



3 THE SIMULATIONS 

The code used for the simulations is a significantly improved 
version of the TreeSPH code we used for our previous work 
on galaxy formation (SGP03). A similar version of the code 
has been used recently to simulate clusters of galaxies, and 
a detailed description can be found in Romeo et al. (2006). 
Here we briefly mention its main features and the upgrades 
over the previous version of SGP03 — see also Sommer- 
Larsen (2006). 

(i) The basic equations are integrated by incorporat- 
ing the "conservative" entropy equation solving scheme of 
Springel & Hernquist (2002), which improves the numerical 
accuracy in lower resolution regions. 

(ii) Cold high-density gas is turned into stars in a prob- 
abilistic way as described in SGP03. In a star-formation 
event an SPH particle is converted fully into a star parti- 
cle. Non-instantaneous recycling of gas and heavy elements 
is described through probabilistic "decay" of star particles 
back to SPH particles as discussed by Lia et al. (2002a). In 
a decay event a star particle is converted fully into a SPH 
particle, so that the number of baryonic particles in the sim- 
ulation is conserved. 

(iii) Non-instantaneous chemical evolution tracing 10 el- 
ements (H, He, C, N, O, Mg, Si, S, Ca and Fe) has been 
incorporated in the code following Lia et al. (2002a, b); the 
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Figure 4. Spatial (projected) distributions of gas-phase oxygen in the "medium-density" field galaxy region, at z=3, for the K98 
simulation. Results for 5 different temperature bins, corresponding to log(T)=4, 4.5, .., 6, are shown. All temperature bins are represented 
using about 5000 "dots", which for each temperature bin correspond to equal amounts of oxygen mass. Hence, each figure illustrates 
directly the spatial distribution of the gas-phase oxygen mass, not gas mass. The cold-phase (log(T)=4) is shown by red dots in all 4 
figures. In addition the log(T)=4.5, 5, 5.5, and 6 phases are shown by blue dots in the corresponding panels. 



algorithm includes supernova of type II and type la, and 
mass loss and chemical enrichment from stars of all masses. 
Most of the simulations presented in this paper have been 
undertaken using three different Initial Mass Functions: the 
Kroupa (1998) IMF (denoted K98 in the following), de- 
rived for field stars in the Solar Neighborhood, the standard 
Salpeter IMF (S), and the more top-heavy Arimoto & Yoshii 
(1987) IMF (AY), which is well suited for the modeling of 
elliptical galaxies, as well as galaxy clusters. More detail is 
given in Lia et al. (2002a, b). 

(iv) Atomic radiative cooling is implemented, depending 
both on the metallicity of the gas (Sutherland & Dopita 
1993) and on the meta-galactic UV field, modeled after 
Haardt & Madau (1996). Moreover, a simplified treatment 
of radiative transfer, switching off the UV field where the gas 
becomes optically thick to Lyman limit photons on scales of 
~ 1 kpc, is invoked. 



(v) Star-burst driven winds are incorporated in the simu- 
lations at early epochs (z>5-6), as strong early feedback is 
crucial to largely overcome the angular momentum problem 
(SGP03). A burst of star formation is modeled in the same 
way as in SGP03: when a star particle is formed, further self- 
propagating star formation is triggered in the surroundings; 
the energy from the resulting, correlated SNII explosions 
is released initially into the interstellar medium as thermal 
energy, and gas cooling is locally halted to reproduce the 
adiabatic super-shell expansion phase; a fraction of the sup- 
plied energy is subsequently converted (by the hydro code 
itself) into kinetic energy of the resulting expanding super- 
winds and/or shells. The super-shell expansion also drives 
the dispersion of the metals produced by type II supernova? 
(while metals produced on longer timescales are restituted 
to the gaseous phase by the "decay" of the corresponding 
star particles, see point 2 above). 
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Figure 5. Spatial (projected) distributions of log(T)=6 gas-phase oxygen in the "medium-density" field galaxy region, at z=3, for the 
K98 (top left), Salpeter (top right), AY (bottom right) simulations. Also shown is the result for a K98 simulation of a proto-elliptical 
region (bottom left) — see Fig. 4 for more detail 



At later epochs, only a fraction (typically, 20%) of the 
stars induce efficient feedback, and star formation is no 
longer self-propagating so that no strong star-bursts are 
triggered by correlated SN explosions. This allows the 
smooth settling of the disc (see SGP03 for all details). Star 
formation efficiencies have been chosen such that realistic 
disc galaxy gas fractions result at z~0 (SGP03), and can 
also be shown to match the Kennicutt (1998) star forma- 
tion relation quite well. 

AGN driven feedback has not been invoked in the simula- 
tions, as it is unlikely to play a major role in the formation 
of, at least, disc galaxies (e.g., Sommer-Larsen 2006) — the 
most common galaxy type in the Universe. The higher stel- 
lar UV escape fractions as well as LyC luminosities found 
by Razoumov & Sommer-Larsen (2007) at z = 3.6, com- 
pared to z < 3 , also support the notion that massive stars 
in galaxies become progressively more important sources of 
ionizing photons as one goes back in time, as the comoving 
number density of quasars declines rapidly at z > 3 (e.g., 



Richards et al. 2006). This provides additional circumstan- 
tial evidence that AGN feedback is of minor importance, at 
least at z > 3, the relevant redshift range in this paper. 

The galaxies were drawn and re-simulated from a 10/i -1 Mpc 
box-length dark matter (DM)-only cosmological simulation, 
based on the "standard" flat A Cold Dark Matter cosmo- 
logical model (h = 0.65, Qo = 0.3, as = 1.0); our choice 
of h and as is slightly different from presently more popu- 
lar values (0.7 and 0.9 respectively), but this has little im- 
pact on the resulting galaxy properties (e.g., Portinari & 
Sommer-Larsen 2007). When re-simulating with the hydro- 
code, baryonic particles were "added" to the original DM 
ones, which were split according to an adopted baryon frac- 
tion /(, = 0.15. The gravity softening lengths were fixed in 
physical coordinates from z=6 to z~0 and in co-moving co- 
ordinates at earlier times. 

The simulations are run with resolutions of msPH = 
m. = 3.6xl0 3 -4.5xl0 7 M Q ,m DM = 2.0x 10 4 -3.3xl0 7 M© 
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Figure 6. Comparison of predicted and actual gas-phase oxygen mass distributions for combinations of various regions simulated using 
the K98 IMF at high resolution. The top left plot shows the result for the combination of three regions, each containing a galaxy simulated 
using the "zoom-in" technique plus a number of companions. The three main galaxies have a mean My of —14.8. The red solid histogram 
shows the actual oxygen mass distribution of the three regions combined. The blue dashed histogram shows what the method outlined 
in the text predicts on the basis of the combined galaxy luminosity function of the three regions. The top right plot shows a similar 
comparison for a combination of 5 high-resolution regions containing main galaxies of an average My of —19.5. The bottom right plot 
is for the combination of the low and medium density "field" regions described in the text, containing a total of 103 resolved galaxies 
of My down to about —22.5. Finally, the bottom left plot shows the results for a higher density "proto-elliptical" region, containing 60 
galaxies of My down to about -22.5. 



and £sph = e. = 0.08 — 3.8 kpc going from the smallest in- 
dividual galaxy simulations to the largest proto-cluster sim- 
ulation. Each region containing an individual proto-galaxy 
was simulated using 0.2-2.2 million particles in total; in 
addition two "low" and "medium" galaxy density regions 
were simulated using 1.4 and 1.6 million particles, a proto- 
elliptical region using 1.4 million particles, and a proto- 
cluster region using 2.3 million particles (the high resolution 
"Virgo" simulation of Romeo et al. 2006). 



Images of some of the simulated galaxies are available 
at http://www.tac.dk/~jslarsen. 



4 THE APPROACH 

The approach adopted is, based on the simulations, to 
characterize a given galaxy at 2=3 and of absolute mag- 
nitude My by a distribution of ISM/IGM oxygen mass 
m (Mv,log(T)), where log(T)=3.5, 4, 4.5 , ..,7, and a given 
value of log(T) corresponds to gas temperatures in the range 

T = [ W lo S (T)-0.2 5AQl o g (T)+ .2 5] K Almogt nQ Qxygen ig 

found at temperatures lower than 10 3 ' 5 K, mainly because 
the radiative cooling function is effectively truncated below 
about 8000 K. If molecule formation and molecular cooling 
had been invoked in the hydro/gravity simulations, part of 
this gas would be in a much colder, high-density, predom- 
inantly H2 phase (see also Sec. 8.1). Moreover, it is found 
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Figure 7. Similar to Fig. 6, but for simulations using the AY IMF. The top left plot shows the result for the combination of three 
regions, each containing a galaxy simulated using the "zoom-in" technique plus a number of companions. The three main galaxies have 
a mean My of —14.0. The top right plot shows results for a combination of 2 high-resolution regions containing main galaxies of an 
average My of —16.3. The bottom right plot is for the combination of the low and medium density "field" regions described in the text, 
containing a total of 40 resolved galaxies of My down to about —21.5. Finally, the bottom left plot shows the results for a very high 
density "proto-cluster" region, containing 86 galaxies of My down to about —25.5. The black dotted histogram shows the prediction using 
a constant f approach (see text); the blue dashed histogram shows the prediction resulting from making the My < —22 f modification 
described in the text 



that the log(T)=3.5 as well as 4 ISM/IGM oxygen is spa- 
tially associated with the stellar galaxies (similar spatial ex- 
tents), so the log(T)=3.5 and 4 bins have been merged, and 
are denoted by log(T)=4. Hence, all gas colder than 10 4 K 
is included in this phase. Moreover, for field galaxies at z~3 
almost no oxygen is found in the log(T)>7 phases, which 
will therefore be ignored in the following. 

Once the functions mo(My, log(T)) have been de- 
termined for a given IMF, the cosmologically averaged 
(gaseous) oxygen temperature distribution can be deter- 
mined by folding these functions with the (dust corrected) 
rest-frame V-band luminosity function at z ~ 3. We stress 
that mo{My, log(T)) pertains only to the galaxy itself, not 
including any of its satellite galaxies — see further below. 



To determine mo{My, log(T)) we proceed as follows: 
in Fig. 3 (left) is shown the cumulative mass of ISM/IGM 
oxygen, Mo(r; My, log(T)), around three (mostly) isolated 
galaxies of My = —16.1, —21.4 and —25.3, respectively, sim- 
ulated using the K98 IMF. Fig. 3 (right) is similar, but for 
the same three galaxies, of My = —13.8, —19.3 and —25.4, 
simulated with the AY IMF For all three galaxies the re- 
gion from r v i T to 8 r„; r is shown 1 , and the different curves 
correspond to log(T)=4, 4.5 , ..,7. As can be seen from the 



In this paper the virial radius, r v ; r , is the radius corresponding 
to, at z=3, an over-density 181 times the mean matter density of 
the Universe at this redshift, as appropriate for a top— hat collapse 
in the adopted ACDM cosmology (e.g., Bryan & Norman 1998). 
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Figure 8. Observational and extinction corrected, modified V- 
band restframe galaxy luminosity functions, based on K-band ob- 
servations by Shapley et al. (2001). Red dashed curve shows the 
original Shapley et al. LF of faint end slope a = —1.85, red solid 
curve the corresponding extinction corrected, modified LF. Blue 
dashed curve shows the result of fitting a LF of faint end slope 
a = —1.57 to the Shapley et al. data, blue solid curve the corre- 
sponding extinction corrected, modified LF. 



Figure 9. Cosmologically averaged gas-phase oxygen density 
temperature distributions for the three IMFs considered. Solid 
and dashed, green histograms correspond to the K98 IMF and 
the extinction corrected, modified galaxy luminosity functions of 
faint end slopes -1.85 and -1.57, respectively. The blue and red 
histograms show the corresponding quantities for the Salpeter 
and AY IMFs, respectively. For comparison, a density of 10 7 
Mq/(1i~ 1 Mpc) -3 corresponds to a density parameter of Ho = 
po/pcTit = 2.5 x 10~ 5 , for h=0.7 



figure, for a given galaxy and log(T), the cumulative gas oxy- 
gen mass gently increases with r by factors ~1.5-4, over the 
radial range shown. The one exception is the log(T)=4 case, 
where the gas metals are mainly spatially correlated with 
the central galaxies, and the curves are quite flat (with the 
exception of the very bright galaxy, which has been selected 
from the proto-cluster simulations). We stress that the mass 
of gas-phase oxygen around a given (isolated) galaxy has 
been contributed not only by the stars in the galaxy itself, 
but also by stars in all its satellite galaxies (and in principle 
also by "inter-galactic" stars, as will be discussed below). 

We select a sample of "base" galaxies consisting of 
(mainly) isolated galaxies ranging in absolute V-band mag- 
nitude from about —12 to — 25. It is assumed that in or- 
der to account for the mo(My, log(T)) of a given base 
galaxy one should include all gas oxygen mass inside of a ra- 
dius f(M v ,log(T)) = r(My,\og(T))/r vir (My), to be deter- 
mined in the following, such that mo(My, \og(T))=Mo(f x 
r V ir', My; log(T)). "Isolated" galaxies are in this paper de- 
fined as galaxies not containing companions of stellar mass 
larger than 1/3 of the mass of the galaxy itself within 8 r v i r . 
The very brightest galaxies, however, of My < — 23, were 
selected from the higher density proto-cluster simulations, 
and the criterion was relaxed to not containing companions 
of stellar mass larger than 1/2 of the mass of the galaxy 
itself within 4 r„; r . The relation between r„; r and absolute 
V-band (rest-frame) magnitude My is shown in Fig. 2 for 
the three IMFs under consideration. 

For a given galaxy, f will be finite for two reasons: a) 
no galaxy is perfectly isolated, such that, even for isolated 
galaxies, with increasing f increasing amounts of gas-phase 
oxygen associated with other (smaller) companion galaxies 



are included, and b) for certain values of log(T) (depending 
on the IMF) galaxies are found to "share" a common pool 
of gas-phase oxygen. In practice it is found that for none of 
the IMFs considered, and any log(T), does f exceed 8. For 
the brightest galaxies, My < — 23, f does not exceed 4. 

Fig. 4 shows, for log(T)=4.5, 5, 5.5 and 6, the spatial 
distribution of the gas-phase oxygen for a region containing 
52 resolved galaxies with My in the range —15 to —22.5, 
simulated using the K98 IMF. Each dot represents, for a 
given log(T'), an equal amount of gas-phase oxygen mass, so 
the figure displays the actual distribution of gas-phase oxy- 
gen mass, not gas mass. As can be seen, the spatial oxygen 
mass distributions depend strongly on the gas temperature. 

In Fig. 5 is shown, for the three IMFs considered, the 
distribution of log(T)=6 oxygen mass. As can be seen, the 
oxygen mass distribution also depends strongly on the IMF 
adopted. For comparison is also shown the corresponding 
plot for a K98 simulation of a proto-elliptical, high-density 
region, which is biased towards many large galaxies (bottom 
left) It is seen how the log(T)=6 oxygen mass distribution 
in this case is concentrated near the large galaxies, more so 
than for the more representative "field" galaxy case. 

We now approximate f (Mv, log(T)) by f(log(T)), i.e., 
independent of My. We show further below that this is a 
good approximation, except for the very brightest galax- 
ies, for which a correction has to be made. For each IMF, 
f(log(T)) is determined as follows: i) for a number of re- 
gions containing up to about 200 resolved galaxies we a) 
determine the absolute V-band magnitudes of the N galax- 
ies in the region, Mv,i, Mv,n, and b) the total gas-phase 
oxygen mass temperature distribution of the region 
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fo , region 



•J regior. 



po(x; log(T)) dx 



(2) 



ii) for a given set of base galaxies, and given values of the 
f(log(T))'s, one can predict the total mass of oxygen as a 
function of log(T) as 



/S"0og(T)) = E£ 1 M^ d (fxr„ lI .(Mv,0;M v , 1 ;log(T)) ,(3) 



where Mq~ c (f x 7w; My,i\ log(T)) is determined by linear 
interpolation in Mv as well as f(log(T)) (f=0,l,...8) using 
the corresponding values for the base galaxies. The values 
of f (log(T)) are then obtained such as to provide the best 
simultaneous match of /g^° n d (log(T)) to /§ < :^ g 1 on (log(r)) 
for a number of simulated regions, using standard least 
squares fitting. 

For all three IMFs this is done for the same two regions, 
representing "low-density" and "medium- density" "field" 
environments, with galaxies of Mv > —22.5. The two regions 
are simulated using a total of 1.4 and 1.6 million particles, 
respectively. For the AY IMF, two additional regions cen- 
tered on smaller galaxies are also used (contrary to the K98 
and S IMF cases, for the AY simulations such small galaxies 
have substantial amounts of not only log(T)=4-6, but also 
log(T)=6.5 and 7 gas-phase oxygen associated) — these re- 
gions are each simulated using about 300,000 particles in 
total. Finally, to test the approach also in a proto-cluster 
environment, for all three IMFs the same such region was 
simulated using a total of 2.3 million particles (the high- 
resolution "Virgo" region of Romeo et al. 2006) . 

As will be shown in the next section, in particular for 
log(T)=4.5, 5 and 5.5, a significant contribution to the cos- 
mic gas-phase oxygen budget comes from small galaxies with 
M v > - 17 to -19. To check that the f (log(T)) fits also work 
in this range, comparison was made to combined data-sets 
of regions with Mv of the central galaxy ranging from —14 
to —19, all simulated at high resolution, with typically 200- 
300 000 particles per region. The results for the K98 IMF are 
shown in Fig. 6, including also the result for a proto-elliptical 
region. In general, the predicted distributions match the ac- 
tual ones very well, except for the very smallest galaxies, 
where the predicted distribution falls somewhat short of the 
actual one in the log(T)=4, 4.5 and 5 temperature bins. 
This implies that the results obtained in the following sec- 
tion for the faintest galaxies are actually lower limits. Given 
the small magnitude of the effect we shall, however, neglect 
this in the following. 

In Fig. 7 a similar comparison is shown for AY IMF 
simulations. Again, the match between predicted and actual 
distributions is, in general, good. For the smallest galaxies 
the distributions are somewhat different (see above). For the 
proto-cluster region, the predicted distribution exceeds the 
actual one, with /o^c g ion d 0°g(^)) being up to about twice 
as large as /ofro|ion(log(T)). This is due to the very large 
galaxy density in such environments. Using this region to 
solve for an alternative set of f (log(T)) results in f(log(T)) 
values of about half the ones obtained as described above. 
As the main contribution to fo^cglon(^°s(T)) is associated 
with the larger galaxies in the region, My ~ —24 to —25.5, 
/o^region(log(T')) in all regions can be matched by a model, 
where f(log(T)) changes linearly from the "field" values to 
the "proto-cluster" values between Mv = —22 to —24, as 
shown in Fig. 6. However, in any case this correction for 



the very large galaxies is of no consequence to the main 
results presented in the following section. This is because the 
contribution from such galaxies to the overall cosmic metal 
budget is very minor, as will be detailed in the following 
section. 

Finally we note, that even in typical "field" environ- 
ments, some fraction of the stellar systems originally formed 
have subsequently been disrupted through tidal stripping 
and other dynamical processes. The gas-phase metals pro- 
duced by these "inter-galactic" stars are "automatically" ac- 
counted for by the above approach, but in any case the frac- 
tion of such stars is quite small, <5-10% (paper I). 



5 THE TEMPERATURE DISTRIBUTION OF 
THE COSMIC GAS-PHASE OXYGEN 

The average cosmic density of gas-phase oxygen at z=3, as 
a function of gas temperature, can now be obtained as 



(po(log(T))> 



dN(L\ 
dLv 



■mo(M v (Lv),log(T))dL v ,(4) 



where dN(Lv)/dLv is the rest-frame V-band, extinction 
corrected galaxy luminosity function at z=3. 

In paper I the "true" z ~3 galaxy luminosity function 
was determined by extinction correcting a modified version 
of the observationally determined luminosity function (LF) 
of Shapley et al. (2001) of faint end slope a = -1.85. The 
original Shapley et al. (2001) LF, as well as the extinction 
corrected, modified LF, are shown in Fig. 8. Also shown are 
the corresponding LFs obtained by fitting a LF of faint end 
slope a = —1.57 (see below) to the Shapley et al. (2001) 
data. The faint end slopes of these z ~ 3 LFs are steeper 
than what is found for z ~ LFs. Results of semi-analytical 
galaxy formation models (e.g., Lacey et al. 2005) indicate 
that this is a natural outcome of hierarchical strucure for- 
mation scenarios. 

Using the extinction corrected, modified Shapley 
et al. LF in eq. 4, and the mo{My(Lv), log(T)) functions 
determined for the three IMFs as described in the previous 
section, results in the distributions shown in Fig. 9. For the 
K98 IMF, the amount of oxygen in the log(T) >4.5 phase is 
about 1.1 times the amount in the log(T)=4 "cold" phase. 
For the S and AY IMFs, the corresponding numbers arc 
about 1.7 and 3.2, respectively. Hence, in particular, for the 
latter two IMFs, the majority of the gas-phase oxygen is 
in the "warm" /"hot" rather than cold phase. Moreover, for 
the AY IMF, the hot to cold/warm gas-phase oxygen mass 
ratio is 1.1, so the bulk of the oxygen is actually in the hot 
log(T) > 5 phase, rather than in the cold/warm phase. These 
are some of the main results of this paper. 

In Fig. 9 is also shown the results of adopting a galaxy 
luminosity function with less steep faint end slope than the 
Shapley et al. one. Following paper I, results are presented 
for a luminosity function of faint end slope a — —1.57, as 
found for the rest-frame UV z ~3 luminosity function by, 
e.g., Adelberger & Steidel 2000). Although the results for 
this LF are qualitatively similar to the ones described above, 
the amounts of warm and hot phase oxygen decrease some- 
what compared to the cold phase ones. The warm/hot to 
cold gas-phase oxygen mass ratios are for this LF 0.8, 1.1 
and 2.1 for the K98, Salpeter and AY IMFs, respectively. 
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Figure 10. Contribution to cosmologically averaged gas-phase oxygen density as a function of galaxy V-band absolute magnitude for 
the K98 (left) and AY (right) IMFs. Curves corresponding to different gas temperatures are colour coded as in Fig. 3 (log(T)=4: green, 
4.5: blue, 5: cyan, 5.5: magenta, 6: black, 6.5: red). The results shown are for the extinction corrected modified Shaplcy ct al. (a = —1.85) 
galaxy luminosity function. For comparison with the cosmological density parameter, see caption of Fig. 9. 




Figure 11. Fractional contribution to the cosmologically averaged gas-phase oxygen density by galaxies brighter than My, for the 
extinction corrected modified Shapley et al. luminosity function (left) and the luminosity function of faint end slope -1.57 (right). Green, 
blue and red curves correspond to results for the K98, Salpeter and AY IMFs, respectively. Long-dashed curves correspond to gas-phase 
oxygen, short-dashed to stellar oxygen, and solid to the total. 



Moreover, for the AY IMF, the hot to cold/warm gas-phase 
oxygen mass ratio is unity. 

An obvious question is now: for the different gas tem- 
peratures, which galaxies contribute mainly to the cosmic 
gas-phase oxygen budget? In Fig. 10 this is shown using the 
modified Shapley et al. luminosity function, for the K98 and 
AY IMFs, respectively. It is seen that for the K98 IMF the 
main contribution to the log(T)=4.5 phase comes from small 
galaxies, My ~ —16, and for the AY IMF this is the case 
for the log(T)=5 and 5.5 phases as well (the reason for the 
large difference between the Mv-distributions of log(T)=4 
and log(T)=4.5, say, gas-phase metals is that, whereas the 



log(T)=4 gas is of fairly high density and associated with all 
galaxies, log(T)=4.5 gas is typically diffuse and associated 
with halos of comparable virial temperature, i.e. fairly small 
galaxies) . 

Given this, and in relation to observational studies, it 
is clearly of interest to determine how large a fraction of the 
gas-phase oxygen is associated with galaxies to a certain lim- 
iting magnitude. In Fig. 11 this is shown by the long-dashed 
curves (for the three IMFs considered) for the a — —1.85 and 
— 1.57 LFs. Typical Lyman Break galaxy metallicity deter- 
minations at 2 ~3 only probe to about one magnitude below 
L* , i.e., My r 22. For the a = -1.85 LF, it is found that 
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less than about 20% of the cosmic gas-phase oxygen will be 
associated with galaxies brighter than this, irrespective of 
the choice of IMF. For the a = —1.57 LF the corresponding 
fraction is about 30%. Hence, most of the cosmic gas-phase 
oxygen is associated with galaxies considerably fainter than 
L* , and is also in this sense "missing" . 

In calculating the results shown in Fig. 11 it has been as- 
sumed that the observational LFs maintain a constant faint 
end slope down to My = —10.5. As can be seen from the 
figure, in particular for the a — —1.85 case this is somewhat 
critical to the shape of the curves shown. At redshifts z~3 
the LBGs luminosity function is only probed down to a few 
magnitudes below L* (but see below). Selecting galaxies us- 
ing the Lya line it is possible to probe 2-3 magnitudes fainter 
still (Fynbo et al. 2001, 2003; Gawiser et al. 2006; Nilsson 
ct al. 2007). This, however, only demonstrates the existence 
of star-forming galaxies at these redshifts with Mv ~ — 17 
to —18, but the data are not good enough to put strong con- 
straints on the shape of the luminosity function. Jakobsson 
ct al. (2005) find that the magnitudes of z > 2 GRB selected 
galaxies are consistent with being drawn from the steep lu- 
minosity function, but the statistics are still poor. We note, 
however, that Bouwens et al. (2007) were actually able to 
probe the UV luminosity function at z ~ 4 — 6 down to 
4-5 magnitudes below L* uv and find steep faint end slopes, 
a 1.7. 

Our calculations show, however, that if the slope of 
the faint end LF flattens significantly at some brighter My, 
then the curves still fairly well represent the result, after ap- 
propriate renormalisation at this limiting Mv (if the faint 
end slope is assumed to be constant to even fainter than 
My = —10.5, the results are not much changed compared to 
the Mv = —10.5 results). Moreover, we note that neither the 
Two-Degree Field Galaxy Redshift Survey (2dFGRS) bj lo- 
cal luminosity function (Norberg et al. 2002), nor the Sloan 
Digital Sky Survey (SDSS) U-band local luminosity func- 
tion (Baldry et al. 2005) show any flattening of the faint end 
slope to 5-6 magnitudes fainter than M*. Moreover, some 
galaxy cluster luminosity functions tend to show a steepen- 
ing faint end slope: Milne et al. (2007) find a steepening of 
the slope in the Coma cluster at a Mr range of about —14 
to —9, and Yamanoi et al. (2007) find a similar result for 
the Hydra I cluster. 

Given the significance of the results obtained in this 
section on the distribution of cosmic gas-phase metals, it is 
important to consider the following two sources of potential 
model dependence. 

First, it is possible that the results depend on the spe- 
cific choice of "base" galaxies. We verified explicitly, how- 
ever, that choices of alternative sets of isolated sets of "base" 
galaxies, and subsequently going through the procedure de- 
scribed in the previous section, lead to very similar results 
as the ones described here. Hence, the results presented are 
very robust to this. 

Second, it would be expected that the amount of (in 
particular) log(T)=4 cold gas and metal mass available in 
a simulation at any given time would depend on the star- 
formation efficiency adopted (Sec. 3). To this end, three K98 
galaxy simulations, of galaxies of z~3 Mv= — 19.1, —21.5, 
and —22.3 and surrounding regions, were run with twice the 
standard star-formation efficiency. Comparing the results of 
these simulations to the corresponding standard ones, it was 



found that the total mass of gas-phase oxygen at z~3 was 
virtually unchanged, and that the mass of log(T)=4 oxygen 
was reduced by <5-10% compared to the standard case. 
Hence, none of the main conclusions about log(T)=4 met- 
als, obtained in this and the following sections, are affected 
by such a change. Moreover, a doubling of the star-formation 
efficiency would lead to too small 2=0 disc galaxy gas frac- 
tions, and likely also too small disc galaxy angular momenta 
(SGP03), likely too small 2^3 DLA cross-sections (Ellison 
et al. 2007, Sommer-Larsen et al. 2007), and star-formation 
rates above the Kennicutt (1998) relation. 



6 C IV PROPERTIES OF THE LOW TO 
INTERMEDIATE DENSITY IGM 

Interesting tests of the abundance and thermal properties 
of the low to intermediate density IGM are provided by 
OSO absorption line studies of the Lyman metal forest, as 
probed by Civ. Songaila (2001) first determined f2(Civ) 
by integrating the total column density of systems between 
10 12 < iV(Clv)< 10 15 cm -2 — see also Schaye et al. (2003). 
As shown by, e.g., OD06, this range of column densities trace 
moderate IGM gas over-densities, 5 = p ga s/p g as ~ 1 — 100, 
or fin ~ 10~ 5 — 10~ 3 cm~ 3 , at 2~3. To obtain a correspond- 
ing estimate of fi(C iv) based on our simulations we adopted 
the following approach: 

To convert from total C (as given in our simulations) 
to C IV abundances, ionisation corrections are required. We 
base our conversion on ionisation corrections determined by 
OD06 using CLOUDY (Ferland et al. 1998). OD06 assumed 
the Haardt & Madau (2001) 2~3 UV background field di- 
vided by a factor 1.6 at all redshifts, in order to match the 
observed mean Lya flux decrement, DL ya - This may seem 
inconsistent with our use of the Haardt & Madau (1996) 
UVB in the cosmological simulations, but as shown by Croft 
et al. (1998), because photo-ionisation is sub-dominant in 
gas dynamics such a correction yields virtually identical re- 
sults as having done the simulations with the above (re- 
duced) background. The resulting ionisation corrections 
were kindly supplied to us by Ben Oppenheimer and Romeel 
Dave, in the form of look-up tables. 

In Fig. 12 (left) is shown the cumulative Q(C iv) as 
a function of 5. For the two simulated "field" galaxy re- 
gions the average ratio of the cumulative C iv mass to to- 
tal (gas-phase) oxygen mass has been determined. This has 
then been multiplied by Qo,gas, determined on the basis 
of the results obtained in the previous section. As is clear 
from this section, the value of !)o,jos depends on the actual 
form of galaxy luminosity function. To obtain lower and up- 
per limits on Qo,gas, we evaluate it at My=— 21 and —10, 
respectively, assuming the a = —1.85 luminosity function. 
My=— 21 corresponds approximately to the detection limit 
of Shapley et al. (2001), whereas Mv = — 10 corresponds to 
assuming that the faint end slope remains constant to this 
very faint magnitude. As can be seen from Fig. 11, using the 
a — —1.57 LF provides less conservative limits. In the fig- 
ure is shown the result of applying the lower and upper lim- 
its for the K98, Salpeter and AY IMFs, respectively. Using 
other simulated regions, centered on smaller galaxies, yields 
similar results. Using the proto-cluster region to determine 
the average ratio of the cumulative C iv mass to total (gas- 
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Figure 12. Cumulative cosmic Civ density, in units of the critical, Q(Clv) as a function of the gas over-density 8 = p ga s//5gas (left). 
Green, blue and red curves correspond to the K98, Salpeter and AY IMFs. Solid curves corresponds to upper limits, dashed to lower. 
Also shown, as a horizontal, black solid line, is the median of ten observational estimates at redshifts in the range 1.5-4.5 taken from 
Songaila (2001), Boksenberg, Sargent & Rauch (2003) and Songaila (2005). Horizontal black dashed lines show 1-<t deviations, where a 
is the variance of the observational estimates. The plot to the right shows the median C IV abundance weighted temperature of the IGM 
as a function of gas over-density, for the three IMFs. 



phase) oxygen mass results in lower values of fi(C iv), due to 
the temperature dependence of the C/C IV ionisation correc- 
tion combined with the higher average IGM temperature of 
the proto-cluster region. However, proto-cluster regions are 
very special sites, and, as can be seen from the results of the 
previous section, the contribution to the cosmic metal pro- 
duction associated with the very luminous galaxies typical of 
such regions, My ~ —25, is very small (in the proto-cluster 
regions containing 100-200 galaxies, almost half of the gas- 
phase metals are associated with five galaxies of My < — 24) . 
Also shown in the figure is the median value of 10 obser- 
vational estimates in the redshift range z=2-4, where the 
run of observed fi(Crv) with redshift is essentially flat (e.g., 
Songaila 2001, Schaye et al. 2003, Songaila 2005), as well as 
the observational variance. The observational estimates were 
taken from Songaila (2001), Boksenberg, Sargent & Rauch 
(2003) and Songaila (2005). 

As can be seen, the chemical yield of the K98 IMF ap- 
pears too low to match the observational estimates, whereas 
the Salpeter and AY IMFs can match the observations, and 
an IMF with a yield in between these two appears optimal. 
We caution, however, that our estimates of fi(C iv) depend 
on the adopted gas-phase oxygen to Civ mass conversions. 
The optimal would be to carry out simulations of > 100 
h~ Mpc box size cosmological volumes at the (high) resolu- 
tion of our galaxy formation runs, but as mentioned previ- 
ously, this is currently computationally prohibitive. Hence, 
our result above on the cosmic IMF should be seen as indica- 
tive only. The main point is that both the Salpeter and AY 
IMF simulations can match the observational constraints on 
fi(Clv). 

Observations of C iv line widths for the above column 
density range can be used to probe the thermal properties 
of the low-density IGM. Three components contribute to 
the line widths: a) thermal broadening, b) spatial broaden- 



ing due to Hubble expansion across the physical extent of 
the absorber, and c) turbulent broadening. The last compo- 
nent can be ignored, since the IGM at these over-densities is 
very quiescent (Rauch et al. 2005). Observational estimates 
of Civ line widths by Boksenberg, Sargent & Rauch (2003) 
indicate b ~10 km/s in the range z ~1.5-4.5 for 10 13 < 
Af(Clv)< 10 14 cm -2 absorbers. This can be used to derive 
upper limits to the C iv abundance weighted IGM temper- 
ature. As the therm al com ponent of the C iv line width is 
given by b th = 3.7^/T/lO 4 , Civ abundance weighted IGM 
temperatures of less than ~5x 10 4 K are indicated. In Fig. 12 
(right) we show C iv abundance weighted IGM temperatures 
a function of S in the "field" galaxy regions for the three 
IMFs considered (results for other regions are quite similar) . 
It is seen that models in general satisfy the above IGM tem- 
perature criterion. Moreover, our predictions of the abun- 
dance weighted IGM temperatures agree well with the pre- 
dictions of OD06 for their best fitting "momentum-driven" 
wind model. 

Finally we note, that in the future it will be possible 
to probe the metal enrichment of the IGM to even lower 
densities using Ovi — see, e.g., Schaye et al. (2000). 



7 COMBINING THE COSMIC STELLAR AND 
GAS-PHASE OXYGEN DISTRIBUTIONS 

With the results on the cumulative cosmic gas-phase oxygen 
distributions presented in the previous section, and the re- 
sults on the similar distributions for the stellar oxygen from 
paper I, we can now assess the relative importance of the two 
components. In Fig. 11 results for the normalised gas-phase 
oxygen distributions are shown by long-dashed curves, the 
stellar distributions are shown by the short-dashed curves, 
and the combined distributions by the solid curves (for the 
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Figure 13. Cumulative fractional contributions of oxygen mass in stars and gas, up to a given My, are shown by the green solid lines, 
separating the red and blue hatched regions (stellar /gas-phase oxygen). All results are for the a = —1.85 LF, and the plots shown are 
for the K98 (top left), Salpeter (top right) and AY (lower right) IMFs. Also shown in these plots, by black solid lines, is the fractional 
partition between oxygen mass in stars and gas at any given My. Finally, in the lower left plot is shown the cumulative cosmic oxygen 
density in gas+stars (solid lines) and stars only (dashed lines). In this plot the IMFs are colour coded as follows: green: K98, blue: 
Salpeter, red: AY. For comparison with the cosmological density parameter, see caption of Fig. 9. 



three IMFs considered) for the a = —1.85 and -1.57 LFs (we 
note that stellar oxygen refers to oxygen in galactic stars 
only, but since the amount of oxygen in intergalactic stars 
is, in comparison, very small (paper I), we neglect this com- 
ponent in the following). 

As can be seen, the combined distributions are for all 
three IMFs dominated by the gas-phase oxygen, and in- 
creasingly so going from the K98 to the AY IMF. For the 
a = —1.85 LF, it is found that less than about 25% of 
the total cosmic oxygen is associated with galaxies brighter 
than My = —22, irrespective of the choice of IMF. For the 
a — —1.57 LF the corresponding fraction is about 35% for 
the AY IMF, and about 40% for the K98 and Salpeter IMFs. 

Fig. 13 shows, for the a — —1.85 luminosity function, 
and the three IMFs, the cumulative partition between gas- 
phase oxygen and stellar oxygen as a function of My. It is 



seen, that if the faint end slope of -1.85 holds down to My ~ 
— 12, then the cumulative gas-phase oxygen fractions are 72, 
79 and 92% for the K98, Salpeter and AY IMFs, respectively, 
so, as stated above, for all three IMFs, the amount gas-phase 
oxygen dominates over the amount stellar oxygen. For the 
a = —1.57 luminosity function these fractions drop slightly 
to 64, 71 and 88%. 

As discussed in Sec. 5 the calculations assume a con- 
stant faint end slope down to Mv=— 10.5. If the luminosity 
function is assumed to display a significant flattening at a 
brighter magnitude than this, then the figure can still be 
used to determine the gas-phase/stellar oxygen partition to 
such a limiting My. Assuming, for example, that the faint 
end slope is constant to at least 6 magnitudes below M* 
(i.e., Mv~ — 17), which is the case for local luminosity func- 
tions (Sec. 5), results in lower limits on the cumulative gas- 
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phase oxygen fractions of 65, 72 and 90% for the a — —1.85 
luminosity function, i.e., not much different from the results 
quoted above. 

We can now finally compare the integral constraints 
on the average cosmic oxygen density at z~3, obtained in 
Sec. 2, to what is obtained by combining detailed high- 
resolution galaxy formation simulations with the observed 
(corrected) z~3, V-band luminosity function. In Fig. 14 
we show the cumulative oxygen density versus My for the 
three IMFs and the two faint end luminosity function slopes 
considered. Also shown are the constraints from the "min- 
imum", "median" and "maximum" oxygen production rate 
density models discussed in Sec. 2. In order to enable a 
meaningful comparison, the latter results should be com- 
pared to the former evaluated at certain limiting magni- 
tudes. These magnitudes should be chosen to be brighter 
than or equal to about My — —21.0, —20.5 and —20.0 for 
the K98, Salpeter and AY IMFs, respectively, for the follow- 
ing reasons: 

Kcnnicutt (1998) finds a relation between UV luminos- 
ity and star formation rate (SFR), viz. 

SFR (Moyr' 1 ) = 1.4 x 10" 28 Amf Mergs s _1 Hz _1 ), (5) 

with /3imf = 1-0 for the Salpeter IMF. Given that the yield 
of the Kroupa (1998) IMF is smaller than for the Salpeter 
IMF, and vice versa for the Arimoto-Yoshii IMF, /3imf will 
be larger and smaller than unity for these two IMFs, respec- 
tively. Quantitatively we find that /3imf —1.7 and 0.4 for the 
K98 and AY IMFs, respectively. 

The absolute UV magnitude (A ~1500 A) characteris- 
ing the observed galaxy luminosity function at z~3 (i.e., 
corresponding to L^j Vz=3 ) is MJ Vz=3 ~ —21.0 (e.g., Stei- 
del et al. 1999, after correction to the adopted cosmology, 
Sawicki & Thomson 2006). Following Sawicki & Thomson 
(2006), and using eq. 5 above, one can show that this abso- 
lute magnitude corresponds to a (un-extincted) star forma- 
tion rate of about 15 /3imf Mq/yt. With a dust correction 
of about a factor 5.5 at z~3 (Sec. 2) the above observed 
My V z=3 hence corresponds to a true (un-obscured) SFR of 
about 82 /3imf Mq/jt. We shall now consider two models for 
the A -1500 A extinction as a function of redshift for z >3. 
Model A assumes a constant extinction factor of 5.5 at all 
z >3 (cf. Sec. 2), whereas model B is a low-extinction model 
assuming factors of 4.2, 3, 2 and 1.5 at z=3, 4, 5 and 6, re- 
spectively (cf. Bouwens et al. 2006). For model A, a galaxy 
of observed Ljjv = 0.1xLyy z=3 will be characterised by a 
SFR~8 /3imf Mq /yr. From our large sample of galaxy mod- 
els we find that galaxies with such star formation rates at 
any z >3 will have an (un-extincted) My at 2=3 brighter 
than about -22.0, -21.5 and -20.5 for the K98, Salpeter 
and AY IMFs, respectively. For model B, the above Myy z=3 
corresponds to a true SFR of about 63 /3imf Mq/yt, and 
galaxies of observed Luv = 0.1xZf£/v,*=3 wn l have SFRs 
of about 6.3, 4.5, 3.0 and 2.3 Mq/jx at z=3, 4, 5 and 6, 
respectively. We find that only galaxies of (un-extincted) 
My at z=3 brighter than about —21.5, —21.0 and —20.5 for 
the above three IMFs, respectively, will satisfy this. More- 
over, even assuming (very conservatively) zero extinction at 
z > 6, corresponding to a "limiting" (un-extincted) SFR of 
about 1.5 Mo/yr, all galaxies of such z>6 SFRs will be be 
brighter than about M v =-21.5, -21.0 and -20.5 at z=3 
for the three IMFs, respectively. Assuming lower luminosity 




Figure 14. Contribution to the cosmologically averaged total 
(gas-phase and stellar) oxygen density by galaxies brighter than 
My, for the modified Shapley et al. luminosity function, and the 
luminosity function of faint end slope -1.57. Green, blue and red 
curves correspond to results for the K98, Salpeter and AY IMFs, 
respectively, (for comparison with the cosmological density pa- 
rameter, see caption of Fig. 9). In the right side of the plot, the 
a = —1.57 curves are above the corresponding a = —1.85 ones. 
Also shown, by horizontal dashed lines, are the constraints from 
the integrated cosmic oxygen production history derived in Sec. 2. 
Finally, the lower V-band luminosity limits, at which the predic- 
tions for the various IMFs and LFs should be compared to the 
integral constraints, are shown by green, blue and red vertical 
dashed lines for the K98, Salpeter and AY IMFs, respectively — 
see text for details. 



limits of Af v =-21.0, -20.5 and -20.0 for the three IMFs 
is hence very conservative — these limits are indicated in 
Fig. 14 by vertical dashed lines. 

From Fig. 14 it follows that galaxy models based on the 
K98 IMF can not meet the constraint set by the observed 
UV luminosity density history — the oxygen (and general 
metal) yield of this IMF is simply not sufficiently large. The 
same is the case for the Salpeter IMF, though the a — —1.57 
model comes close to matching the lower bound of 0.32 x 10 r 
M©/(/i -1 Mpc) 3 (see also below). On the other hand, mod- 
els based on the Arimoto-Yoshii IMF, match the integral 
constraint well. 

Summarising, one of the main results obtained in this 
paper is that galaxy formation models based on the Kroupa 
(1998) IMF (and any other IMF of similar chemical yield) 
are strongly excluded by the cosmic enrichment history con- 
straint. Other arguments, why the average cosmic IMF must 
have a larger chemical yield than what is typical for a so- 
lar neighborhood one, have been given by, e.g., Portinari 
et al. (2004), D'Antona & Caloi (2004), Serjeant & Harrison 
(2005), Lucatello et al. (2005), Loewenstein (2006), Prant- 
zos & Charbonnel (2006) and Weidner & Kroupa (2006), 
but see also Elmegreen (2006). 

Perhaps even more interestingly, models based on the 
Salpeter IMF, the arguably most extensively used model of 
the stellar IMF, are also excluded, though the a = —1.57 
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models only marginally so. Given that the Salpeter IMF is 
routinely used in translating UV luminosities into star for- 
mation rates for high redshift galaxies, this is also an impor- 
tant result. Although it would be inappropriate to assign a 
high statistical significance to this result, given all the uncer- 
tainties involved, it is obvious from Fig. 14 that the models 
based on the more top-heavy Arimoto-Yoshii IMF provide 
a much better match to the median value for the integral 
constraint. This is in particular the case, since the above 
lower absolute luminosity limits of My=- 21.0, —20.5 and 
—20.0 for the three IMFs, respectively, are likely to be very 
conservative. 



8 MISSING METALS AND COMPARISON TO 
OTHER WORKS 

In Sec. 5 it was shown that for the a = —1.85 lumi- 
nosity function, less than about 1/4 of the cosmic oxygen 
is associated with Lyman Break galaxies sufficiently bright, 
My > — 22, for direct abundance determination, using oxy- 
gen lines emitted from HII regions around young stars in 
the galaxies (e.g., Pettini et al. 2001). For the a = —1.57 
luminosity function this fraction increases slightly, to about 
35%. Furthermore, in Sec. 7 it was shown that for both lu- 
minosity functions the major part of the cosmic oxygen is 
in the gas-phase, rather than in stars. In particular, for the 
Salpeter and Arimoto-Yoshii IMFs, which emerge from the 
previous section as the more plausible, the gas-phase oxygen 
fraction exceeds 70%. 

FSB05 found a factor of about five discrepancy between 
the amount of metal in the stars of Lyman Break galaxies 
and what is predicted from the integrated UV luminosity 
density history. They assumed typical Lyman Break galaxy 
stellar masses of 2 x 10 10 M Q , corresponding to Mv~ —22.5, 
i.e., close to My z=3 . 

In Fig. 13 is shown, for the a = —1.85 luminosity func- 
tion, the cumulative stellar cosmic oxygen density. For the 
K98 IMF, the ratio between the stellar oxygen density to 
Mv = —21 and what is obtained from the median model 
(Sec. 7) is 0.10. For the Salpeter and AY IMFs, evaluated 
at My = —20.5 and —20, respectively, the corresponding 
ratios are 0.17 and 0.20. For the a = —1.57 luminosity func- 
tion the corresponding ratios are 0.14, 0.24 and 0.29, for the 
K98, Salpeter and AY IMFs, respectively. Given that the 
AY models over-predict the metallicities of Lyman Break 
galaxies of My ~ —22 to —23 by about 0.2 dex relative to 
observations (paper I), the above ratios for the AY models 
should be reduced to 0.13 and 0.18 for the a — —1.85 and 
-1.57 luminosity functions, respectively. 

The above results are obtained, however, by including 
stellar oxygen mass all the way down to My~ —20.5. If one 
only includes stellar oxygen mass to Mv~ —22.5, all the 
above fractions are reduced by about a factor of two. Hence 
the discrepancy discussed above, denoted by FSB05 and oth- 
ers as "the missing metals problem" , is actually about twice 
larger than originally found by FSB05. 

DO06 predicted that at z~3, about 50% of the gas- 
phase metals should reside in the diffuse IGM (gas outside 
the virial radii of galaxy halos, and of T < 3xl0 4 K). We 
can not compare the results found in the previous sections 
directly to theirs, since different gas-phase criteria have been 



used. However, our log(T)=4.5 gas criterion is similar to 
theirs for diffuse IGM, though we also include log(T)=4.5 
metals inside of galaxy virial radii in our estimate (see also 
below; note also that we find almost all the log(T)=4 phase 
metals to reside in gas of fairly high density, tih > 0.1 cm~ 3 , 
so there is essentially no overlap between this phase and 
DO06's diffuse phase). 

For the a = -1.85 LF we find log(T)=4.5 to total gas 
phase metal fractions of 34, 31 and 26% for the K98, Salpeter 
and AY IMFs, respectively. For the a — —1.57 LF, the cor- 
responding fractions are 21, 20 and 18%. If we try to mimick 
the criterion of DO06 better, by selecting all metals in gas 
of T<3xl0 4 K and <5 = p gas / p gas <W0, then for the field 
galaxy regions we find "diffuse IGM" oxygen fractions of 
4, 4 and 6% for the K98, Salpeter and AY IMFs, respec- 
tively. If we use the proto-cluster regions, which are the 
only simulations at our disposal with a resolution compa- 
rable to the (fairly modest) numerical resolution of DO06's 
32 /i _1 Mpc box size simulation, the corresponding fractions 
drop to ~1% for all IMFs. However, we stress that the tem- 
perature of the proto-cluster IGM is larger than that of the 
average IGM, but in any case significantly lower "diffuse 
IGM" metal fractions, than predicted by DO06, are indi- 
cated. 

DO06 used the Springel & Hernquist (2003) "sub-grid" 
approach to model star-formation. They predict that about 
20% of the gas-phase metals reside in star-forming gas. In 
our models, which invoke explicit two-phase modeling of the 
ISM, gas of T~10 4 K and n H > 0.1 cm~ 3 is potentially star- 
forming. Hence we compare our results for the log(T)=4 
gas-phase metal content to the above of DO06 (as noted 
above, almost all log(T)=4 metals reside in gas of uh ^0.1 
cm -3 ). For the a = -1.85 LF we find log(T)=4 to total 
gas phase metal fractions of 47, 37 and 26% for the K98, 
Salpeter and AY IMFs, respectively. For the a = —1.57 LF, 
the corresponding fractions are 57, 47 and 32%. Hence we 
find somewhat larger metal fractions in star-forming gas, 
than do DO06. A more important difference is, however, 
that DO06 at z~3 find no metal containing gas at densities 
uh cm -3 . This result is strongly at variance with our 

results, and it would seem difficult for such models to match 
the large amount of DLA cross section observationally de- 
tected at such redshifts, although further analysis obviously 
is required to clarify this. 

Finally, we can compare the z~3 gas phase metal tem- 
perature distribution predicted by DO06 to our results. 
Qualitatively, our predictions for the K98 and Salpeter IMFs 
agree with that of DO06, in yielding temperature distribu- 
tions, which peak at T~10 4 K and decrease towards larger 
temperatures. However, for the AY IMF the results are very 
different. Moreover, quantitatively our results for the K98 
and Salpeter IMFs disagree significantly with those of DO06 
at the "high-T" end: DO06 find that < 5% of the gas-phase 
metals reside in gas of T>10 4 ' 75 K. In comparison, we find 
for the a = —1.85 LF that the corresponding fractions are 
20, 32 and 51% for the K98, Salpeter and AY IMFs, respec- 
tively, and for the a = -1.57 LF, 21, 32 and 50%. Hence for 
any of the IMFs considered, we predict significantly larger 
fractions of "high-T" gas-phase metals, than do DO06. 
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Figure 15. Left: Fractional cumulative oxygen mass in the absorbers above log(N(HI))=20.3 versus the cumulative absorber area, 
starting at the highest HI surface densities. Red and blue curves correspond to the proto-galaxies of 2=0 V c =245 km/s and 180 km/s, 
respectively. "Face-on" projections are shown by solid and "edge-on" by dashed curves. Center: Same, but versus HI column density. 
Right: Oxygen and iron abundance versus HI column density; for clarity only "face-on" projections are shown. Red and blue curves show 
[O/H], green and cyan [Fe/H]. 



8.1 Metals in the cold gas-phase and DLA 
absorbers 

The log(T)=4 phase is found to be the most prominent 
metal containing gas-phase for the K98 and Salpeter IMFs, 
and also to be quite prominent for the AY IMF. As the met- 
als in this phase typically reside in gas of densities hh j^O.l 
cm~ 3 , one would expect a significant fraction of these met- 
als to be situated at column densities typical of DLAs, viz. 
N(HI)>10 20,3 cm -2 . However, in general DLAs have very 
low abundances - significantly less than Galactic stars with 
similar ages (Pettini et al. 1990, 1994; Lu et al. 1996; Kulka- 
rni & Fall 2002; Prochaska et al. 2003; Akerman et al. 2005; 
Zwaan et al. 2005; Erni et al. 2006). Expressed in terms of 
the global metal content of the LBG stars, we predict similar 
to three times as large amounts of metals in the log(T)=4 
phase (depending on the IMF), whereas observations of 
DLAs only indicate a fraction of about 10-20% (e.g., FSB05, 
Prochaska et al. 2006). It is obviously of importance to un- 
derstand the reason for this apparent discrepancy. Although 
we will defer a thorough discussion of DLAs to a forthcom- 
ing paper (Sommer-Larsen et al. 2007), we briefly address 
the above issue in the following. 

Ellison et al. (2007) analysed two of the very high res- 
olution simulations, described in the following section, in 
relation to DLA properties. The two simulations represent 
the formation and evolution of two disc galaxies, of z=0 
characteristic circular velocities 14=245 and 180 km/s. The 
galaxies were selected from a larger sample to represent two 
different disc formation evolutionary paths: for the 14=180 
km/s galaxy the disc starts growing steadily already by 
z~2.5, whereas for the other galaxy disc growth is merger 
induced, and the disc grows strongly from z~l to (see 
also Robertson et al. 2004). Ellison et al. determined the 
DLA/sub-DLA characteristics of these two proto-galaxies 
at 2=3.6, 3.0 and 2.3, with focus on determining neutral 
column density distributions, and the probability of detect- 
ing coincident 100 kpc-scale DLA/sub-DLA absorption in 
individual galaxies at such redshifts — we refer the reader 



to Ellison et al. for more detail, as well as images of the 
objects. 

Here we build on the analysis by Ellison et al. of these 
two (proto-)galaxies. In Fig. 15 (left) we show the nor- 
malised cumulative oxygen mass in the absorber above 
log(N(HI))=20.3 versus the cumulative absorber area, start- 
ing at the highest HI surface densities, log(N(HI))~24, and 
going down to log(N(HI))=20.3. For each galaxy results 
for "face-on" and one "edge-on" projections are shown. 
As can be seen, about 90% of the total oxygen mass of 
the absorber resides in 1% of the total absorber area for 
both projections. From Fig. 15 (center) it is moreover seen 
that this oxygen mass is associated with column densities 
log(N(HI)) > 22.5. DLAs of such high column densities have 
never been detected in QSO spectra, which of course is 
not surprising given the above findings. Moreover, at such 
high densities, formation of molecular hydrogen is likely to 
take place (Schaye 2001b; Zwaan & Prochaska 2006). Al- 
though the process of H2 formation is not included in the hy- 
dro/gravity simulations, it is clear that the properties of the 
high density gas can be significantly affected by H2 forma- 
tion (e.g., Pelupessy et al. 2006). Greve & Sommer-Larsen 
(2007) showed that the effects of H2 formation can be ap- 
proximately determined post-process on the basis of the cos- 
mological galaxy formation simulations — this will be one of 
the main topics of a forthcoming paper on DLAs (Sommer- 
Larsen et al. 2007). H2 molecules have been detected in 
DLAs, but only at relatively low fractions of less than a tenth 
relative to atomic hydrogen (Ledoux et al. 2003). Moreover, 
H2 is preferentially detected in the highest metallicity DLAs 
(Petitjean et al. 2006). At this point it is sufficient to note 
that a) only a very small fraction of QSO sight-lines will 
probe the high- metallicity regions of DLAs, and b) due to H2 
formation these regions may possibly be difficult to probe in 
neutral hydrogen (Zwaan & Prochaska 2006). We also note 
from Fig. 15 (right) that the metal abundances of the high 
column density regions are significantly larger than what is 
typically observed in DLAs, whereas at log(N(HI)) < 22 they 
are not (e.g., Pettini et al. 2003). Finally, we note that the 
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high column density regions will likely be characterised by 
significant dust contents, which could further bias against 
selecting QSO sight-lines passing through such regions. Ev- 
idence for dust in DLAs has been reported (Pei, Fall & 
Bechtold 1991; Vladilo et al. 2006). This may naturally ex- 
plain why DLA column densities of log(N(HI)) > 22 are very 
rare (Vladilo & Peroux 2005, Vladilo et al. 2006). However, 
studies of radio selected DLAs (free from dust-bias) have 
found similar column density and metallicity distributions 
as for optically selected DLA samples (Ellison et al. 2001, 
2004; Ellison, Hall & Lira 2005; Akerman et al. 2005; Jor- 
genson et al. 2006). Also, Murphy & Liske (2004) find no 
excess absorption towards QSOs with DLAs (but this study 
is based on an optically selected DLA sample). It remains 
to be clarified to which extent dust bias is important for 
DLA column density and metallicity studies, especially at 
log(N(HI)) > 22. In particular, a larger sample of radio se- 
lected DLAs is needed (Jorgenson et al. 2006, Ellison 2007, 
private communication) . 

We note that DLAs at z ~ 2 with near solar metallicity 
have been found (Ledoux et al. 2002). Also, the Sloan Digi- 
tal Sky Survey has been used to strongly increase the sample 
of DLAs (Prochaska et al. 2005) and among these 5% have 
very strong metal- line absorption (Herbert-Fort et al. 2006). 
Some of these 5% are metal rich DLAs, although some could 
be relatively metal poor DLAs with very large hydrogen col- 
umn densities. About 60% of the metal-strong DLAs have 
detected molecular hydrogen (Milutinovic et al. in prepara- 
tion). 

Interestingly, DLAs with log(N(HI)) > 22 have been de- 
tected in the spectra of the optical afterglows of Gamma- 
Ray Bursts (Watson et al. 2006; Jakobsson et al. 2006; 
Prochaska et al. 2007). The selection function for GRB 
sight-lines is completely different than the cross-section se- 
lection of QSO DLAs. GRB sight-lines are strongly cor- 
related with the blue light of their host galaxies (Bloom, 
Kulkarni & Djorgovski 2002; Fruchter et al. 2006) most 
likely tracing the location of the most massive stars 
(M > 2OM0 , Larsson et al. 2007) . Significant dust obscura- 
tion is observed for at least some of the log(N(HI)) > 22 GRB 
DLAs (Watson et al. 2006). Furthermore, it has been found 
that the metallicities inferred for DLAs in GRB hosts are 
systematically higher than for QSO DLAs and that GRBs 
hence may offer a new probe of early cosmic chemical evo- 
lution complementary to the QSO DLAs and LBGs (Fynbo 
et al. 2006; Savaglio 2006; Prochaska et al. 2007). 

In conclusion, the low metallicities inferred from QSO 
DLAs do not exclude that significant amounts of metals can 
reside in the log(T) = 4 phase. The reason for this is most 
likely that the total cross-section for this component is only 
of order a few percent of the total DLA cross-section (see also 
Johansson & Efstathiou 2006) and that the current DLA 
samlpes are too small to uncover them (see also Zwaan & 
Prochaska 2006), but dust bias, as well as H2 formation, will 
probably also be of importance. 



9 NUMERICAL RESOLUTION 

In SGP03 and Sommer-Larsen (2006) it has been shown 
that the results of cosmological galaxy formation and evo- 
lution using our code are, in general, robust to changes of 
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Figure 16. Gas-phase oxygen temperature distributions for three 
proto-galactic regions simulated at both very high (blue dashed 
histograms) and normal (red solid histograms) numerical reso- 
lution. The three proto-galaxies become at 2=0 disc galaxies of 
characteristic circular velocities V c =245, 180 and 66 km/s, re- 
spectively. The oxygen masses for the smaller galaxy have been 
multiplied by a factor 50 for clarity. The slight offset between 
the results for the smaller galaxy is due to high resolution proto- 
galactic region being slightly smaller than the corresponding nor- 
mal resolution one. 



the numerical resolution. However, in this paper we present 
gas-phase metal temperature distributions and other results, 
which have not been presented before. It is hence clearly 
important to demonstrate that these results are resolution 
robust as well. To this end we carried out three very high 
resolution K98 IMF simulations, of eight times higher mass 
and two times higher force resolution compared to "stan- 
dard". We simulated three (proto-)disc galaxies, which at 
2=0 have characteristic circular velocities of 245, 180 and 
66 km/s. The galaxy forming regions were represented by a 
total of about 2.2, 1.2 and 1.3 million particles, compared 
to the about 290 000, 150 000 and 210 000 particles used at 
"normal" resolution (for computational reasons the region 
simulated at very high resolution for the last galaxy was 
slightly smaller than the corresponding "normal" resolution 
one). Gas and star particle masses were m gas =m*=9.1x 10 4 
and 2.5xl0 3 /i _1 Mq for the first two galaxies and the last 
one, respectively. The corresponding dark matter particle 
masses were moM=5.2x 10 5 and 1.4xl0 4 /i _1 M Q . More- 
over, gravitational (spline) softening lengths of e gas =e* = 190 
and £dm=340 ft~ 1 pc, respectively, were adopted for the 
first two galaxies, and 57 and 100 fe~ 1 pc for the latter. 
Other results for the first two runs have already been pre- 
sented in Sommer-Larsen (2006), Razoumov & Sommer- 
Larsen (2006), Greve & Sommer-Larsen (2007), Laursen & 
Sommer-Larsen (2007) and Ellison et al. (2007). 

In Fig. 16 is shown the gas-phase oxygen temperature 
distributions for the three very high resolution runs, together 
with the results of the normal resolution ones — note that 
the oxygen masses for the smaller galaxy have been multi- 
plied by a factor 50 for clarity. It is seen that the general 
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agreement is very good; the slight disagreement for the small 
galaxy is due to the fact that the region simulated at very 
high resolution was slightly smaller than the corresponding 
normal resolution one. 

Due to computational limitations we did not perform 
similar resolution tests for the two other IMFs considered, 
but we have no reason to believe that simulations based on 
these IMFs would perform less well in resolution tests. We 
conclude that the results obtained in this paper are robust 
to resolution changes. 



10 SUMMARY AND CONCLUSIONS 

The global temperature distribution of the cosmic gas-phase 
oxygen at z~3 has been determined by combining high reso- 
lution cosmological simulations of individual proto-galactic, 
as well as larger, regions with extinction-corrected, observa- 
tionally based, V-band (rest-frame) galaxy luminosity func- 
tions (LFs) of faint end slopes a = —1.85 and —1.57. 
The simulations have been performed with three different 
stellar initial mass functions (IMFs), a Kroupa (K98), a 
Salpeter (S) and an Arimoto-Yoshii (AY), spanning a range 
of a factor of five in chemical yield and specific SNII en- 
ergy feedback. Gas-phase oxygen is binned according to 
T as log(T)<4.0 ("cold"), log(T) ~ 4.5 ("warm"), and 
log(T) ~5.0, 5.5, 6.0, 6.5, 7.0 ("hot" phases). Below we sum- 
marize results for the a = —1.85 LF, but results for the 
a = —1.57 LF are similar. 

Oxygen is found to be distributed over all T phases, 
in particular for the ("top-heavy") AY IMF. But, at vari- 
ance with previous works, it is found that, for the K98 and 
Salpeter IMFs, the most important phase is the cold one, 
which contains 47 and 37% of all gas-phase oxygen, mainly 
in gas at fairly high density, riH >0.1 cm~ 3 , and potentially 
star-forming. Moreover, the cold phase alone contains 1.3, 
1.5 and 3.2 times the mass of oxygen in galactic stars for 
the three IMFs. The implications of this in relation to ob- 
servational damped Lyman-a absorber (DLA) studies are 
discussed on the basis of very-high resolution simulations of 
two (proto-)disc galaxies, with emphasis on oxygen and iron 
abundances. It is concluded, that the reason why current 
DLA surveys only detect a cold ISM metal fraction of about 
20% relative to the metal mass in galactic stars is that the 
total cross-section for the high-metallicity component is only 
of order a few percent of the total DLA cross-section, and 
that the current DLA samlpes are too small to uncover it. 
Moreover, in addition, dust bias, as well as H2 formation, 
will likely also be of importance. 

In relation to "missing metals" it is found that the ratio 
of gas-phase to stellar oxygen mass is 2.7, 3.9 and 13, and 
the ratio of warm+hot to cold gas-phase oxygen mass is 1.1, 
1.7 and 3.2 for the three IMFs. For the AY IMF, the hot 
phases actually contain more oxygen than the cold+warm. 

In conclusion, a significant fraction of the cosmic oxy- 
gen may be difficult or impossible to detect. In addition, 
it is found that less than about 20-30% of the cosmic oxy- 
gen will be associated with galaxies brighter than Mv ~ 
—22, i.e., the faintest galaxy luminosities probed by current 
LBG metallicity determinations (about one mag. below L* ). 
Hence, 70-80% of the cosmic oxygen is also in this sense 
"missing" . 



From the LBG based, A~1500 A UV luminosity density 
history at z>3, we obtain an essentially IMF independent 
constraint on the mean cosmic oxygen density at 2=3. We 
compare this to what is obtained from our models, for the 
three different IMFs. We find that the (solar neighbourhood 
type) K98 IMF is strongly excluded, as the chemical yield is 
simply too small, the Salpeter is marginally excluded, and 
the AY matches the constraint well. The optimal IMF would 
have a yield intermediate between the S and AY. The K98 
IMF can only match the data if the A~1500 A extinction 
corrections have been overestimated by factor of ~4, which 
seems highly unlikely, cf. Rcddy & Steidel (2004). 

Using carbon abundances, and C to C IV ionisation cor- 
rections, we estimate 0(C iv) at moderate IGM gas over- 
densities, 5 = Pgas/pgas ~ 1 — 100 for the three IMFs, and 
compare to observational results. As above, we find that the 
yield of the K98 IMF is too small to match the data, whereas 
models based on the Salpeter and AY IMFs can match the 
data, with the optimal IMF in between the two. Moreover, 
we show that for all IMFs, C IV abundance weighted IGM 
temperatures are moderate, T<4xl0 4 K, consistent with 
observational constraints on C IV line widths. 
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